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ABSTRACT 


We  are  investigating  Pn  propagation  in  Western  China  based  on  regional  events  recorded  by  the  Hi-CLIMB  (An 
Integrated  Study  of  the  Himalayan-Tibetan  Continental  Lithosphere  during  Mountain  Building)  array  in  the  region. 
Seismic  attributes,  including  arrival  times,  Hilbert  envelope  amplitudes  and  instantaneous  and  speetral  frequencies, 
are  being  used  to  investigate  how  the  erustal  and  upper  mantle  velocity  and  attenuation  structure  affects  the 
propagation  of  Pn  arrivals  in  Tibet.  We  have  eonstrueted  more  than  30  high-quality  regional  seismie  profiles,  and  of 
these,  14  events  have  been  selected  with  excellent  erustal  and  Pn  arrivals  for  further  analysis.  Regional  travel-times 
recorded  by  the  Hi-CLIMB  array  are  used  to  constrain  the  larger  seale  velocity  structure  in  the  region,  with  four 
events  near  the  array  used  to  constrain  the  crustal  structure.  We  are  also  independently  constraining  the  variation  of 
upper  mantle  velocity  and  depth  of  the  Moho  by  modeling  SF-eoupled  shear  waves  at  tclcseismic  distances.  The 
amplitude  of  the  SsPmP  phase  is  particularly  sensitive  to  upper  mantle  velocities  when  it  is  near  critical.  The  initial 
modeling  of  the  SsPtnP  phase  indicates  that  the  Moho  beneath  southern  Lhasa  is  over  75  km  deep  with  Pn  velocities 
greater  than  8  km/s.  In  contrast,  the  data  sampling  the  Qiangtang  terrane  north  of  the  Bangong-Nujiang  suture  shows 
thinner  crust  with  Pn  velocities  less  than  8  km/s.  This  is  generally  consistent  with  the  results  from  the  Pn  travel- 
times  for  regional  events  recorded  on  the  Hi-CLIMB  array.  Nonetheless,  there  is  additional  variability  of  the 
regional  travel-time  data,  suggesting  further  3D  complexities.  Seismic  amplitude  and  frequency  attributes  arc  being 
extracted  from  the  erustal  and  Pn  wave-trains,  and  these  data  are  being  compared  with  numerical  and  analytical 
results  for  models  with  upper  mantle  velocity  gradients,  which  can  strongly  affect  Pn  amplitudes  and  frequencies. 
The  numerical  modeling  is  performed  using  approximate  ray  and  beam  methods  and  the  more  complete  speetral 
element  method  (SEM).  The  results  from  the  SEM  method  are  in  good  agreement  with  analytical  and  reflectivity 
results  for  different  models  with  upper  mantle  velocity  gradients.  Using  the  SEM  method,  we  are  testing  the  effeets 
of  random  erustal  and  Moho  depth  heterogeneities  to  see  how  this  variability  impacts  the  range  of  Pn  attributes 
observed.  Ultimately  the  approach  we  are  follow  ing  is  a  mostly  data-driven  analysis  of  Pn  seismie  attributes  in  order 
to  investigate  the  trade-offs  of  large  seale  lithospheric  structure,  small  seale  heterogeneities  and  attenuation  in  Tibet. 
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OBJECTIVES 

We  are  analyzing  the  propagation  of  regional  Pn  waves  recorded  by  Hi-CLIMB  in  Tibet.  Considering  the 
station-spaeing  (3  to  8  km)  and  large  aperture  (800  km)  of  the  Hi-CLIMB  seismic  experiment,  our  objective  is  to  use 
these  densely  recorded  data  to  eonstruet  a  self-eonsistent  model  of  Pn  propagation  over  a  large  region  in  western 
China.  In  turn,  a  better  understanding  of  the  effeets  of  propagation  from  analysis  of  the  data  will  improve  frequency 
dependent  attenuation  models  of  regional  Pn  phases,  leading  to  a  elearer  separation  of  elastic  and  anelastie  effeets  in 
a  laterally  varying  erust  and  upper  mantle,  eventually  eliminating  the  assumption  of  a  simplified,  frequeney- 
independent  geometric  spreading  component. 

To  account  for  the  intricate  nature  of  Pn  propagation,  we  are  applying  fast  ray  and  beam  methods  (Nowack  and 
Staey,  2002),  as  well  as  SEM  (Komatitseh  and  Vilotte,  1998;  Komatitsch  et  al.,  2005).  In  contrast,  alternatives  such 
as  the  reflectivity  method,  are  restricted  to  1 D  media.  SEM  ean  be  applied  in  2D  or  3D  and  implemented  in  parallel 
using  MPI  on  large  cluster  computing.  Ray  and  beam  methods  are  faster  and  ean  provide  asymptotic  results 
appropriate  when  a  large  number  of  forward  calculations  are  required,  such  as  in  large-seale  inversions.  However, 
for  small-scale  scattering,  complete  numerical  methods,  such  as  SEM,  may  be  more  appropriate.  As  a  practical 
measure,  essential  features  of  Pn  w  aveforms  are  being  distilled  into  seismie  attributes  ineluding  arrival  times, 
envelopes  of  wave  amplitudes,  and  pulse  frequencies  for  modeling  and  inversion. 

RESEARCH  ACCOMPLISHED 

We  have  constructed  more  than  30  high-quality  seismie  profiles  recorded  on  the  Tibet  Hi-CLIMB  array  and  from 
these  14,  have  been  selected  for  further  analysis.  Similarities  in  the  seismic  profiles  for  events  in  dusters  are 
outstanding.  About  two  thirds  of  the  larger  earthquakes  have  fault  plane  solutions  available  in  the  global  Harvard 
Centroid  Moment  Tensor  (CMT)  eatalog  (Dziewonski  et  al.,  1981)  and  the  report  by  Baur  (2007),  For  those  without 
fault  plane  solutions,  we  have  identified  the  ones  that  ean  potentially  be  resolved  using  regional  waveforms  recorded 
by  Hi-CLIMB  stations  for  their  appropriate  azimuthal  coverage,  or  by  using  teleseismic  waveforms  (Figure  1 ).  The 
analysis  of  moment  tensor  inversion  is  in  progress.  The  fault  plane  solutions  that  eannot  be  obtained  through  this 
method  may  be  roughly  estimated  based  on  geologie  structures  or  historic  solutions. 

As  a  complimentary  approach  to  the  analysis  of  regional  seismograms,  we  are  independently  constraining  the 
variation  of  upper  mantle  velocity  and  the  depth  of  Moho  by  modeling  shear  waveforms  at  teleseismic  distances 
shorter  than  50  degrees.  In  essence,  the  amplitude  of  the  SF-coupled  P  phase  SsPmp  is  very  sensitive  to  the  upper 
mantle  velocities  when  the  rays  hit  the  Moho  near  the  critical  angle.  Our  results  indicate  that  the  significant  decrease 
in  amplitude  of  SsPmp  near  northern  Hi-CLIMB  stations  ean  be  explained  by  slow  velocity  anomalies  in  the 
uppermost  mantle  beneath  northern  eentral  Tibet. 

Regional  events  recorded  by  the  Tibet  Hi-CLIMB  array  with  locations  approximately  in-line  with  the  array  have 
been  ehosen  for  initial  analysis.  However,  other  events  in  speeifie  terranes  have  also  been  selected  for  later 
processing.  Four  of  the  in-line  events  have  epicenters  less  than  200  km  from  the  Hi-CLIMB  array;  these  events  were 
relocated  with  data  from  the  Hi-CLIMB  array  using  the  earthquake  location  program  HYPQ1NVERSE-2000 
(Klein,  2002).  The  locations  of  selected  in-line  events  are  shown  in  Figure  2,  using  the  EHB  (Engdahl  et  al.,  1998) 
or  the  Preliminary  Determination  of  Epicenters  (PDE)  catalogs  for  all  events  except  for  the  four  nearby  events  (L2a, 
L2d,  L3a  and  L3b),  whieh  were  relocated  using  Hi-CLIMB  data.  These  relocations  improved  the  fit  to  the  travel- 
times,  but  moved  the  epicenters  by  less  than  5  km.  For  more  distant  events  in  Tibet,  a  focal  depth  normalization 
procedure  was  used  sinee  events  on  Tibetan  plateau  are  dominantly  less  than  1  5  km  in  focal  depth  (e.g  ,  Molnar  and 
Chen,  1983;  Chen  and  Molnar,  1983;  Langin  et  al.,  2003). 

The  observed  arrival  times  for  the  four  nearby  events  were  then  used  to  refine  the  erustal  velocity  model  beneath  the 
array  with  the  3D  ray  traeing  eode  CRT  (ver.  6.0)  using  the  approach  described  by  Cerveny  et  al.  (1988),  The 
ealeulated  travel-times  were  compared  with  the  observed  arrival  times  in  order  to  iteratively  refine  the  erustal 
velocity  model.  The  initial  velocity  model  was  based  on  seismie  refraction  analysis  of  controlled  shots  for  the 
INDEPTH  III  experiment  by  Meissner  et  al.  (2004).  It  was  found  that  the  relocated  nearby  events  combined  with 
the  refined  erustal  velocity  model  provide  a  reasonable  fit  to  the  observed  erustal  travel-times.  An  example  of  the 
observed  and  ealeulated  travel-times  for  nearby  event  L3a  is  shown  in  Figure  3. 
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Figure  1.  Map  of  the  Hi-CLIMB  array  with  selected  labeled  events  for  which  seismic  profiles  have  been 

constructed.  The  red  focal  mechanisms  were  obtained  from  the  global  CMT  catalog  and  the  gray 
ones  from  Baur  (2007). 

In  order  to  model  the  Pn  phase  arrival  times,  work  is  currently  being  performed  to  pick  the  Pn  branches  and  refine 
the  Moho  depths.  An  example  of  the  ray  tracing  for  the  Pn  diving  ray  for  the  southern  event  Lla  is  shown  in  Figure 
4.  For  events  such  as  Lla,  located  along  the  Himalayan  front,  the  Pn  arrival  times  are  complicated.  There  exist  rapid 
variations  of  Moho  depth  w  ith  a  shallowing  of  the  Moho  along  the  Himalayan  front  and  south  to  the  Indian  shield. 
The  Lhasa  terrane  in  the  central  Tibetan  plateau  has  a  much  greater  Moho  depth,  and  north  of  the  Bangong-Nuijiang 
suture  in  the  Qiangtang  terrane  of  Tibet,  the  Moho  moderately  shallows.  After  modeling  all  the  Pn  arrival  times,  the 
other  seismic  attributes  will  be  modeled,  including  seismic  amplitudes  and  pulse  frequencies 

Selected  In-line  Regional  Events  with  111-CLIMB  Tibet  Array 


Figure  2.  Selected  regional  events  approximately  in-line  with  the  Hi-CLIMB  Tibet  seismic  array. 

For  examples  of  shot  profiles  with  longer  Pn  branches,  we  have  modeled  the  travel-times  for  event  Lie  to  the  south 
and  event  E3b  to  the  north  of  the  Hi-CLIMB  array.  The  left  plot  of  Figure  5  shows  the  travel-times  for  event 
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Preliminary  Travel-Time  Modeling  for  Event  L3A 
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Figure  3.  Initial  travel-time  modeling  of  Event  L3a  to  constrain  crustal  structure  beneath  the  Hi-CLIMB 
array  in  Tibet. 

Lie,  where  the  circles  are  the  observed  times  and  the  crosses  are  calculated  times  using  the  3D  ray  tracer  CRT  for 
the  preliminary  model  shown  in  the  right  plot.  Similar  to  the  modeling  of  travel-times  for  the  nearby  ev  ents  to  the 
array,  the  crustal  phases  for  distances  <260  km  for  event  Lie  are  relatively  slow  down  to  the  Moho.  The  modeling 
of  the  Pn  branch  shows  a  deepening  of  the  Moho  from  the  Himalayan  front  dow  n  to  75  km  for  the  Lhasa  terranc  in 
central  Tibet  with  Pn  speeds  between  8.2  and  8.3  km/s.  Then,  to  the  north  of  the  Bangong-Nujiang  suture  (BNS)  in 
the  Qiantang  terranc,  the  Moho  shallows  by  about  10  km  to  63-65  km  w  ith  slower  Pn  speeds  of  about  7.8  km/s.  In 
the  left  plot  of  Figure  5,  the  Pn  branch  is  first,  resulting  from  the  deepening  Moho,  and  then  becomes  earlier 
resulting  from  the  thinning  of  the  crust  to  the  north.  The  upper  mantle  gradient  inferred  from  the  modeling  is  about 
0.003  to  0.004  km/s  per  km  consistent  with  the  work  of  Philips  et  al.  (2007)  for  the  region. 
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Figure  4.  An  example  of  3D  ray  tracing  of  the  diving  Pn  phase  for  event  Lla  using  a  Tibet  model  with  a 
highly  variable  Moho  topography 

The  second  event-gather  to  be  modeled  is  from  Event  E3b  located  to  the  north  of  the  Hi-CLIMB  array  (shown  in 
Figure  2).  The  travel-times  are  shown  in  the  left  plot  of  Figure  6  where  the  circles  are  the  observed  times  and  the 
crosses  are  the  calculated  times  using  the  3D  ray  tracer  CRT  and  the  preliminary  model  is  shown  in  the  right  plot. 
Again  the  crust  is  relatively  slow  down  to  the  Moho  with  no  indication  of  a  higher  velocity  lower  crust.  The  Pn 
branch  then  arrives  initially  early,  indicating  the  shallower  Moho  to  the  north  and  then  is  delayed,  resulting  from  the 
deeper  Moho  in  the  Lhasa  terrane  in  central  Tibet.  However,  to  match  the  complete  travel-time  curve,  a  slower  Pn 
speed  is  required  to  the  north  with  the  Pn  speed  increasing  to  the  south  as  the  Moho  deepens.  The  modeling  of  this 
event  located  north  of  the  array  is  consistent  with  the  modeling  of  event  Lie  to  the  south  of  the  array.  For  the 
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modeling  of  both  event-gathers,  upper  mantle  gradients  of  between  0.003  and  0.004  km/s  per  km  are  consistent  w  ith 
the  data  and  similar  to  the  results  found  by  Philips  et  al.  (2007)  for  the  region. 


Preliminary  Travel-Time  Modeling  for  Event  LI  E 
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Figure  5.  Initial  travel-time  modeling  of  Event  Lie  to  constrain  upper  mantle  structure  near  the  Hi-CLIMB 
array  in  Tibet.  The  left  plot  shows  the  observed  (circles)  and  calculated  (erosses)  travel-times.  The 
right  figure  shows  the  initial  erustal  and  upper  mantle  model  in  the  vieinity  of  the  Hi-CLIMB  array 
from  travel-time  modeling  for  event  Lie.  Velocities  are  shown  in  km/s  and  single  digits  show  upper 
mantle  velocity  gradients  in  units  of  0.001  (km/s  per  km). 
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Figure  6.  Initial  travel-time  modeling  of  Event  E3b  to  constrain  upper  mantle  structure  near  the  Hi-CLIMB 
array  in  Tibet.  The  left  plot  shows  the  observed  (circles)  and  calculated  (erosses)  travel-times.  The 
right  figure  shows  the  initial  erustal  and  upper  mantle  model  in  the  vieinity  of  the  Hi-CLIMB  array 
from  travel-time  modeling  for  event  Lie.  Velocities  are  shown  in  km/s  and  single  digits  show  upper 
mantle  velocity  gradients  in  units  of  0.001  (km/s  per  km  in  km/s  and  single  digits  show  upper  mantle 
velocity  gradients  in  units  of  .001  (km/s  per  km). 


As  a  preliminary  example  of  seismic  amplitudes,  Figure  7  shows  the  Hilbert  envelope  amplitudes  for  Event  Lie. 

The  dots  in  the  figure  are  envelope  amplitudes  obtained  from  a  two  second  window  after  the  arrival  times.The  lines 
are  curves  proportional  to  1/R  and  1/(R~)  set  to  be  equal  at  200  km.  Between  distances  of  about  260-340  km,  the  Pn 
amplitudes  are  lower  and  then  increase  to  betw  een  340  to  425  km.  After  425  km  and  north  of  the  BNS  in  the 
Qiantang  terrane,  the  Pn  amplitudes  decrease  relative  to  the  two  decay  curves.  Since  the  Moho  decreases  in  depth  by 
about  10  km  across  the  BNS,  the  amplitudes  for  event  Lie  suggest  that  large  scale  structural  features  are  influencing 
the  amplitude  attributes.  Nonetheless,  the  Moho  structure  beneath  the  Himalayas  and  Tibet  is  complicated  and  more 
analysis  is  needed  to  interpret  the  observed  seismic  attributes  on  the  different  seismic  profiles.  We  are  still  analyzing 
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the  results  for  pulse  frequency  attributes,  as  well  as  the  decay  of  amplitudes  with  distance.  However,  we  are 
applying  a  data  driven  approach  to  interpret  the  trade-offs  between  structural  effects  of  the  medium  as  opposed  to 
other  attenuation  effects,  and  not  driven  by  an  assumed  model  for  the  analysis  of  the  observed  seismic  attributes.  In 
any  case,  large-scale  structural  features  need  to  be  taken  into  account  for  the  modeling  of  seismic  attributes,  such  as 
amplitudes  and  pulse  frequencies,  in  addition  to  the  modeling  of  scattering  from  small  scale  features  and  intrinsic 
attenuation. 


First  Arrival  Hilbert  Envelope  Amplitude  for  Event  LIE 


Figure  7.  Hilbert  envelope  amplitudes  of  first  arrivals  for  Event  LI E.  The  two  lines  illustrate  R-l  and  R-2 

decay  curves  shown  to  coincide  with  the  data  at  200  km.  Note  that  after  250  km  the  amplitudes  first 
decay  and  then  get  larger  from  350  to  430  km.  For  distances  greater  than  450  km,  the  amplitudes 
are  lower  again,  possibly  correlated  with  the  Bangong-Nujiang  Suture  in  Tibet. 

We  are  currently  performing  numerical  simulations  of  elastic  wave  solutions  for  later  comparison  with  the  observ  ed 
data.  We  have  implemented  SEM  using  the  parallel  code  Specfem2d  for  the  simulation  of  2D  viscoelastic  and 
anisotropic  media  (Komatitsch  and  Vilottc,  1998;  Komatitsch  ct  al.,  2005).  This  code  allows  us  to  use  the  cluster 
computing  facilities  available  at  Purdue  University  and  the  University  of  Illinois.  These  SEM  calculations  let  us 
compare  with  asymptotic  ray  and  beam  calculations  and  gain  insight  into  the  parameters  that  affect  the  Pn  arrivals, 
including  crustal  and  upper  mantle  velocities  and  gradients,  lateral  variability  in  the  Moho  depth,  and  small  scale 
scattering  and  intrinsic  attenuation  in  the  crust  and  upper  mantle. 


Synthetic  Seismograms  for  Model  1  (Vertical  Component) 


Figure  8a.  Comparison  between  the  results  from  SEM  (dotted)  and  the  reflectivity  results  (solid)  of  Braile 
and  Smith  (1975)  for  a  constant  velocity  30  km  erust  over  a  constant  velocity  mantle.  The 
amplitudes  of  the  reflectivity  results  are  multiplied  by  distance  and  the  SEM  results  are  multiplied 
by  the  square-root  of  distance. 

Figure  8a  shows  an  example  synthetic  record  section  for  a  model  with  a  6.4  km/s  crust,  a  sharp  discontinuity  at  the 
Moho  and  a  constant  8.0  km/s  mantle.  The  times  are  reduced  by  8.0  km/s.  The  PmP  reflected  waves  and  Pn  head 
waves  are  designated  on  the  record  section.  The  solid  lines  are  the  seismograms  computed  using  the  reflectivity 
method  by  Braile  and  Smith  (1975)  and  the  dotted  lines  are  our  calculations  using  SEM  modeling,  where  the 
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amplitudes  are  sealed  by  distance  for  the  reflectivity  calculation  and  by  the  square  root  of  the  distance  for  the  2D 
SEM  synthetics.  The  comparison  between  them  is  very  good,  and  suggests  that  we  are  implementing  the  SEM  code 
correctly.  Note  the  small  amplitude  of  the  pure  head  wave  phase  on  the  synthetics  in  this  ideal  ease. 


Synthetic  Seismograms  for  Model  5  (Vertical  Component) 


Figure  8b.  This  shows  a  comparison  between  the  results  from  our  spectral  clement  calculation  (dotted)  and 
the  reflectivity  results  (solid)  of  Braile  and  Smith  (1975)  for  a  constant  velocity  30  km  crust  over  a 
mantle  with  a  velocity  gradient  in  the  top  5  km.  The  amplitudes  of  the  reflectivity  results  arc 
multiplied  by  distance  and  the  SEM  results  are  multiplied  by  the  square-root  of  distance. 


Figure  8b  shows  an  example  of  synthetic  seismograms  for  a  velocity  model  with  an  upper  mantle  gradient  in  the  first 
5  km  below  the  Moho.  This  model  is  referred  to  as  Moho-5  in  Braile  and  Smith  (1975).  The  solid  seismograms  are 
from  the  reflectivity  calculations  of  Braile  and  Smith  (1975)  and  the  dotted  seismograms  are  from  our  calculations 
using  the  SEM  modeling.  The  reflectivity  amplitudes  are  multiplied  by  distance  and  the  SEM  results  are  multiplied 
by  the  square  root  of  distance.  Again,  the  comparison  betw  een  the  two  results  is  very  good,  suggesting  that  we  are 
implementing  the  SEM  modeling  correctly.  Note  the  large  Pn  amplitudes  resulting  from  the  upper  mantle  gradient  in 
contrast  to  the  ideal  head  wave  ease  in  Figure  8a.  Even  in  the  1970s,  it  was  recognized  that  there  was  a  trade-off 
between  large-scale  structural  features,  such  as  upper  mantle  velocity  gradients  and  other  forms  of  seismic 
attenuation. 
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Figure  9.  The  left  plot  shows  a  comparison  of  the  spectral  element  results  (red  crosses)  with  the  asymptotic 
results  of  Ccrvcny  and  Ravindra  (1971)  for  a  layer  over  a  halfspace  (131  curve)  and  for  a  C+ 
interference  head  wave  for  a  layer  over  a  velocity  gradient.  The  background  curves  are  from 
Ccrvcny  and  Ravindra  (1971).  The  right  plot  shows  the  amplitude  spectra,  where  the  black  curves 
are  from  Ccrvcny  and  Ravindra  (1971)  and  the  red  curves  are  from  the  SEM  results.  Pulse 
distances  for  the  spectra  arc  shown  to  the  right. 
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The  left  plot  of  Figure  9  shows  a  comparison  betw  een  the  log-amplitudes  from  the  SEM  calculations  and  the 
asymptotic  results  arc  from  Cerveny  and  Ravindra  (1971)  for  a  layer  over  a  half  space  and  a  layer  over  a  velocity 
gradient.  The  Moho  is  at  a  depth  of  30  km  with  a  contrast  of  6.4  km/s  over  8  km/s.  For  the  velocity  gradient  case, 
the  gradient  is  0.02264  km/s  per  km.  This  is  larger  than  typically  observed  mantle  gradient  values  (0.001  and  .005 
km/s  per  km)  but  decreases  the  distance  range  for  different  frequency  and  amplitude  effects  to  occur,  making  it  more 
efficient  for  numerical  calculations.  Note  that  the  flattening  of  a  spherical  earth  will  result  in  an  effective  velocity 
gradient  of  about  0.001  km/s  per  km  in  the  flattened  Earth  (Yang  et  al.,  2007).  The  background  figure  is  from 
Cerveny  and  Ravindra  (1971)  where  the  C+  wave  is  the  interference  head  wave,  the  Co  wave  is  the  pure  diving 
wave  and  the  131  wave  is  the  pure  head  wave  case.  The  red  crosses  are  from  the  SEM  calculations  for  the 
interference  and  pure  head  wave  cases.  For  the  asymptotic  results  of  Cerveny  and  Ravindra  (1971),  the  amplitude 
blows  up  at  the  critical  distance  near  80  km  and  the  comparisons  are  shown  only  for  distances  greater  than  1 10  km. 
For  the  C+  interference  head  wave,  the  amplitude  initially  decreases  after  the  critical  distance  and  then  recovers  as 
the  diving  wave  becomes  the  dominant  w  ide-angle  wave.  There  are  some  undulations  in  the  asymptotic  results  of 
Cerveny  and  Ravindra  (1971)  for  distances  greater  than  220  km  for  the  C+  wave  that  are  not  seen  in  the  SEM 
calculations,  but  otherwise  the  results  of  the  SEM  calculations  are  reasonably  consistent. 

The  right  plot  of  Figure  9  shows  the  amplitude  spectra  as  a  function  of  distance  for  the  interference  head  wave.  The 
black  curves  are  from  Cerveny  and  Ravindra  (1971),  and  the  red  curves  are  from  the  SEM  calculations  for  the  same 
layer  over  a  velocity  gradient  model.  Both  results  show  that  the  dominant  frequency  initially  decreases  at  140  km 
w  ith  the  maximum  of  the  amplitude  spectrum  less  than  that  of  the  incident  wave  because  of  the  frequency  decrease 
of  the  pure  head  wave.  As  the  diving  wave  begins  to  emerge,  the  dominant  frequency  increases.  However,  because 
of  tuning  effects  between  the  different  interference  waves,  the  dominant  frequency  overshoots  that  of  the  incident 
wave  to  higher  frequencies.  This  was  also  found  by  Nowack  and  Stacy  (2002)  using  both  Gaussian  beams  and 
reflectivity  calculations. 
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Figure  10.  The  left  Figure  shows  a  comparison  of  the  SEM  results  (red  crosses)  with  the  asymptotic  results  of 
Cerveny  and  Ravindra  (1971)  for  a  layer  over  a  velocity  gradient.  The  red  and  the  blue  crosses 
are  for  the  random  perturbations  on  the  Moho  shown  in  Fig.  11.  The  background  curves  are  from 
Cerveny  and  Ravindra  (1971).  The  right  figure  shows  the  amplitude  spectra  for  the  two 
perturbed  Moho  interfaces  for  a  layer  over  a  velocity  gradient.  The  black  curves  are  from 
Cerveny  and  Ravindra  (1971)  and  the  red  and  blue  curves  are  from  the  SEM  results.  Pulse 
distances  for  the  spectra  are  shown  to  the  right. 

The  left  plot  of  Figure  10  shows  the  amplitude  comparison  between  SEM  calculations  for  two  random  Moho 
interfaces  with  a  standard  deviation  of  0.1  km  and  a  smoothing  length  of  20  km  with  the  unperturbed,  asymptotic 
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results  of  Cerveny  and  Ravindra  ( 1971 )  for  a  layer  over  a  velocity  gradient.  Normal  random  numbers  are  first 
generated  in  Matlab  and  then  are  smoothed  using  three  passes  of  a  boxear  smoothing  of  length  6.5,  6.5  and  7  km  in 
order  to  smooth  out  the  higher  wavenumbers.  The  right  plot  of  Figure  10  shows  the  effect  of  two  random 
perturbations  of  the  Moho  depth  on  the  pulse  spectra  at  different  distances.  The  results  show  that  the  main  features 
of  the  C+  interference  wave  are  preserved  for  small  perturbations  of  the  Moho  depth.  Numerical  results  for  longer 
smoothing  lengths  for  the  perturbation,  but  the  same  average  vertical  deviation,  show  similar  results  but  with  longer 
average  wavenumbers  for  the  amplitude  perturbation.  Numerical  examples  with  a  standard  deviation  of  0.25  km  in 
vertical  depth  perturbation  of  the  Moho  shows  larger  perturbations  of  the  interference  head  wave  amplitude  but  still 
vary  about  the  main  feature  trend  of  the  reference  ease.  Nonetheless,  preliminary  calculations  with  larger  scale 
structure  features  have  more  significant  effeets  on  the  amplitude  and  frequency  eharaeteristies  of  the  Pn  wave  and 
will  be  investigated  in  the  future. 

CONCLUSIONS  AND  RECOMMENDATIONS 

We  have  assembled  a  number  of  record  sections  from  regional  events  recorded  by  the  Hi-CLIMB  array.  For  several 
events  nearby  the  Hi-CLIMB  array,  we  have  relocated  the  events  using  Hi-CLIMB  data.  For  more  distant  events  in 
Tibet,  a  focal  depth  normalization  procedure  has  been  used  sinee  events  on  theTibetan  plateau  are  dominantly  less 
than  15  km  in  focal  depth  (e.g.,  Molnar  and  Chen,  1983).  From  these  record  sections,  seismic  attributes  of  the  Pn 
phase  are  being  extracted  from  the  data,  ineluding  arrival  times,  envelope  amplitudes  and  instantaneous  and  spectral 
pulse  frequencies.  Several  nearby  events  to  the  Hi-CLIMB  array  have  been  used  to  constrain  crustal  velocities  from 
crustal  branches.  For  the  Pn  phase,  the  travel-times  are  being  modeled  using  3D  ray  tracing  in  order  to  constrain 
large  seale  structural  features  including  Moho  depths,  upper  mantle  velocities  and  velocity  gradients.  The  modeling 
of  travel-times  for  all  the  selected  event  gathers  will  result  in  an  integrated  3D  model  for  the  region.  Also, 
teleseismie  analysis  of  the  SsPmp  phase  is  being  performed  to  further  constrain  the  erustal  and  upper  mantle 
structure. 

Numerical  modeling  is  performed  using  the  spectral  element  method  (SEM)  in  order  to  compare  the  effeets  of 
large-scale  erustal  structure  on  seismie  attributes.  For  the  layer  over  a  velocity  gradient  ease,  the  SEM  results 
compare  favorably  with  the  asymptotic  results  of  Cerv  eny  and  Ravindra  ( 1971 )  and  the  reflectivity  calculations  of 
Braile  and  Smith  (1975).  Seismie  amplitude  and  pulse  frequency  results  using  SEM  are  stable  for  small  interface 
perturbations.  Lay  et  al.  (2008)  reported  more  drastic  variations  obtained  from  finite  difference  modeling  for  upper 
mantle  velocity  perturbations.  Presumably  these  results  oeeur  sinee  the  mantle  velocity  perturbations  change  the 
upper  mantle  gradients,  which  eontrol  the  primary  properties  of  the  interference  head  wave.  Nonetheless,  larger 
variations  of  Moho  depth  have  more  pronounced  foeusing  and  defoeusing  effeets  on  Pn  amplitudes  than  shown  here. 
We  are  starting  to  perform  SEM  tests  for  large  seale  structural  models  that  oeeur  in  Tibet  and  these  are  anticipated  to 
have  strong  effeets  on  the  properties  of  the  Pn  wave. 
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